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ABSTRACT 

Assuming that Dark Matter is dominated by WIMPs, it accretes by gravitational attraction and 
scattering over baryonic material and annihilates inside celestial objects, giving rise to a "Dark Lumi- 
nosity" which may potentially affect the evolution of stars. We estimate the Dark Luminosity achieved 
by different kinds of stars in a halo with DM properties characteristic of the ones where the first star 
formation episode occurs. We find that either massive, metal-free and small, galactic-like stars can 
achieve Dark Luminosities comparable or exceeding their nuclear ones. This might have dramatic 
effects over the evolution of the very first stars, known as Population III. 
Subject headings: Population III — early Universe — stars: evolution — dark matter 
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1. INTRODUCTION 

Observations continue to provide evidence for a dark 
component of matter in the Universe, whose motiva- 
tion s and possible ca ndidates are thouroughly reviewed 
by (jBertone et al '05l ). It has been realized quite early 
that scattering between Dark Matter (DM) and baryons 
would result in capture of DM particles over celestial 
bodies, with the most remarkable effect for stars to 
host within their bosom an additional source of en- 
ergy due to the captured DM annihilation. Since then, 
several authors have studied the possibility to exploit 
the effects of an additional, potentially uncxhaustable, 
energy source within stars in order to constraint DM 
prop erties: among the latest , (Moskale nko fc Wai '07f ) 
and dBertone &; Fairbairn '07D study com pact objects; 
(|Fairbairn et al '07l ) and ([Scott et al '07ft have per- 
formed preliminary numerical calculations of the stellar 
evolution. All the results seem to converge over the fact 
that within a neutralino dominated DM scenario, and 
with an astrophysical environment with typical galac- 
tic DM velocity profiles of the order of v «10 2 km/s, 
DM densities required in order to have appreciable effects 
over stellar evolution are p x >10 7 GeV cm -3 , achievable 
only within the central two parsecs of our Galaxy, with 
similar restrictions holding for other types of galaxies: 
DM burning would not affect most of the stars in our 
local Universe. However, so far no one has studied the 
effects of DM burning over the first stars in the Universe: 
according to the current ACDM cosmological description 
the early Universe is quite a peculiar place for stellar for- 
mation to take p lace. The almost total absence of metals, 
(|Iocco et al. '07ft . the small mass and angular momen- 
tum of the young and dense haloes where the first stars 
form, the lack of strong magnetic fields, all intertwine 
together giving the Population III peculiar properties; 
above all, is to be mentioned the charact eristic formation 
of a single, massiv e (3OM -3OOM Q ), (|Abel et al. '02f ) 
and (|Gao et al. '07ft . Massive, metal free stars would 
evolve quickly through the Main Sequence (MS), burn- 
ing all of their fuel within timcscalcs of 10 6 years and 
proceeding quickly through the following clement burn- 
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ing, until at the oxygen burning stage temperature and 
entropy are so high that copious e + e~ production takes 
place, pressure drops and a runaway collapse is triggered. 
This leads to explosive nucleosynthesis which results in 
direct collapse to a Black Hole, or to a direct dramatic 
explosion, or to stellar pulsations follo wed by a final vio- 
lent explosion, (|Heger fc Wooslev '02ft . An "alteration" 
of the basic stel lar formation sce nario which has been 
pointed out by (jSpolvar et al 'Oct ): at times much ear- 
lier than the actual formation of a protostar (T<10 4 K, 
n <10 18 cm~ 3 ) Dark Matter embedded in the halo anni- 
hilates at a rate such that the energy deposited exceeds 
the one loss by the molecular cooling of the gas. Al- 
though these semianalytical results have to be confirmed 
by more detailed numerical simulations it is clear that 
this process could change the history of primordial stel- 
lar formation by preventing the cloud from collapsing, or 
by slowing it down or else inducing fragmentation of the 
baryonic cloud. Although this process is different from 
the one we study it is clearly related to it: when the cloud 
becomes opaque enough to efficiently scatter WIMPs off 
baryons another energy source would be added in the 
center, as opposed to the pure DM annihilation, which 
diffuse everywhere in the halo. Here we present prelimi- 
nary estimates of the DM accretion and annihilation in- 
side a star, as it would result is it had formed in "clas- 
sical" way, without any effect from the DM annihilation 
at lower densities and temperatures. 

2. DARK MATTER CAPTURE AND ANNIHILATION IN 
STARS 

Models of DM capture and annihilation on celestial 
bodies were developed in the 80s by many authors, and 
here we collect their results of most concern for this Let- 
ter. At the equilibrium between ca pture and annihilation 
the event rate is, after ([Gould '87ft : 
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p x is the ambient WIMP energy density, A n is the atomic 
number of the star's nuclei, M n is the nucleus mass, v 
is the WIMP velocity dispersion, is the star's velocity 
with respect to an observer and // = m x /M n , /x_ = Qu — 
l)/2, u(r) is the escape velocity at a given radius r inside 
a star, and subscript * refers to stellar quantities. The 
capture-annihilation process will reach the equilibrium 
on timescales t x with most of the particles concentrated 
within a radius r x , which read respectively: 
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T c and p c are the stellar core temperature an d density, 
respectively, as shown in (|Griest fe Seckel '87t ). Accord- 
ing to them, inside the star the captured DM follows a 
density profile n x , with a central maximum density n£, 
respectively described by: 
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Also, DM particles can transport energy by scatter- 
ing off nuclei inside the star, thus constituting an 
alternative transpor t method inside the star; after 
(|Spergel fc Press '85t ). the quantity of energy trans- 
ported by the WIMPs per unit volume per unit time, 
e x , reads: 
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where the subscript n refers to values for the nucleons. 
L x —m x C depends neither on (av) nor on the ncutralino 
mass m x , whereas they both affect the distribution of the 
DM inside the star and the transport effects. It is worth 
noticing that although the product of an annihilation (of 
the dark matter accreted inside the star), L x depends on 
environmental quantities of dark matter (environment 
density, velocity dispersion) with the typical properties 
of a scattering. Also, both L x and t x are degenerate in 
the product V=<7q x p x /v (actually, plus an additional 
dependence on v, which is a posteriori negligible for this 
illustrative purpose), which we will use as a parameter 
in the following as it conveniently summarizes the DM 
values in the halo, with ctq. 

3. APPLICATION TO EARLY STARS: PRELIMINARY 
RESULTS 

The DM density and its velocity dispersion play a fun- 
damental role in determining the accretion rate over any 



stellar object; in the local Universe, high DM densities 
are achieved only in the very centre of galaxies; in par- 
ticular in the Milky Way, a cuspy profile is believed to 
have formed by accretion around a central Black Hole, 
over the long lifetime o f our galaxy and its ce ntral ob- 
ject, see for instance in (jBertone fc Merritt '05l ) and ref- 
erences therein. DM particles adopt Keplerian velocities 
around the central object, thus reaching their highest ve- 
locities where the DM density peaks; on the other hand, 
most of the stellar mass is not located in the center of 
the Galaxy, but in regions where the DM density is too 
small to consistently affect stellar evolution. 

In a primordial environment the scenario is quite 
different: in a young halo such as the ones where PopIII 
stars form average densities are quite high, as the 
concentration of Dark Matter, c, is higher in a younger 
Universe. Moreover, no central object exists yet so 
DM particles do not have time to develop a keplerian 
velocity profile: the star itself is believed to be the very 
first object to form in the center of the halo (where 
the DM density is highest). In a primordial halo, for 
the whole lifetime of a PopIII star, we can assume 
DM dispersion velocities to be a Maxwell-Boltzmann 
distribution around the central value corresponding 
to the virial temperature of the halo itself. In order 
to estrapolate the values of the DM density p x at the 
location of a primordial star, we quote state of the art 
simulations of the collapse of the the first object. They 
achieve impressive resolutions, having to be stopped at 
central gas (baryon) densities of 0(1O 25 particles/cm 3 ), 
achieved at a radii of «10 9 cm, ([l|); unfortunately this is 
not (yet) enough to resolve the actual star, simulations 
stop following DM particles even earlier, when the 
density is homogenous within a radius of «10 cm, 
thus not allowing us to have more detailed insights on 
smaller scales. For a central region of radius 10 15 cm one 
finds an enclosed mass Mom ~10M q , corresponding to 
a p x «10 12 GeV cm _3 ;this d ata, quote d from the state 
of the art simulations {]]) and (|Turk '08t ). is in agreement 
with the results extrapolated from previous simulations 
and also with the ones obtained by much simpler semi- 
analiti cal models of adi abatically contracted DM halos, 
as in (S polvar et al '08l ). Namely, the halo followed in 
Turk's cosmological simulations has Mfl Q ; o =0(lO 6 M Q ) 
a virial temperature T«10 3 K and therefore a virial 
velocity v «10 6 cm s _1 , according to mass-temperature 
relations evolution with redshift in (jLevine et al. '021 ) 
and references therein. It is likely that the final p x 
in the center of the halo will be higher at the time of 
stellar formation; moreover, the effe cts triggered by a 
process such as the one discussed in (jSpolvar et al '08l ) 
are not trivial to predict. Here we limit ourselves 
to present results using p x =10 12 GeV cm -3 as our 
fiduci al value, adopt ing (av)=3x 10~ 26 cm 3 s _1 af- 
ter (jBertone et al '05h and er =10~ 38 (/10~ 43 )cm 2 
for the spin-dependent(/inde pendent) case, 
(|CDMS '06D / (lSuper-Kamiokande '041 1: this cor- 
responds to a fiducial value for our parameter 
X>=l(T 32 (/l(r 37 )GeV s/cm 2 for the spin depen- 
dent/independent) case. Results for different values 
can be easily rescaled. 

Using the equations introduced in Section [2j we 
have estimated those quantities for WIMPs captured by 
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TABLE 1 

Values for a 75 Mq, initial metallicity Z=10~ 4 , in a 

JEUTRALINO CASE WITH m x = 100GEV. Z>=10~ 32 (/10~ 37 )GeV 
S/CM 2 FOR THE SPIN DEPENDENT(/lNDEPENDENT) CASE. LAST 
LINE REFERS TO THE H SHELL DURING HELIUM BURNING. 
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TABLE 3 

Same as in Table [T] for ZAMS metal-free stars. 
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TABLE 2 

Same as in Table [T] for MS low-mass stars. 



a 75M , extremely metal-p oor (Z=1O~ 4 Z ) star, de- 
scribed in (jWoosley et al.'02h : our results do not change 
within the order of magnitude if we ad opt the values of a 
metal- free, 100M Q star described in (|Marigo et al '01| ). 
In fact, stellar radii of O (R* «1OR ) during the MS 
and R* «lO 3 R cm during the helium burning phases are 
achieved both by the 75M Q star i n (IWooslev et al.'02f) 
and the metal-free IOOMq one in (|Marigo et al 'Oil ), as 
it can be seen from a detailed table in the first case 
and extrapolated from the diagram on page 5 in the sec- 
ond. While this paper was being edited for publication, I 
have received values of Zero Age Main Sequence, metal- 
free stars from P. Marigo (unpublished results after 
(jMarigo et al 'Oil )): Table[3]contains the relevant quanti- 
ties calculated as from equations in Section[2]for these ac- 
tual metal-free stars; r v «10 6 s, nf\=10 19 GeV/cm 3 , and 



x a, "- x —±v GeV/cm 3 
e x «10 6 erg/s/cm 3 for the whole mass range; Lf AMS 
is the nuclear luminosity of the star, i?„ is its radius 



in cm, which we present for direct comparison with 
For the case of the 75M star, in order to estimate the 
mass of the core for the stages following the H burn- 
ing, and therefore the central mass of the star which has 
a definitely differe nt composition, we us e the approxi- 
mate relation from l)Heger fe Wooslev '021 ) (which has to 
be taken with care in this case, as we are out of the 
PISNe range) M ffe =13/24(M ZAM s-20), thus obtaining 
an helium core for the 75M star of Mjj e «3OM . Dur- 
ing the helium burning stage the star's central density 
is p^ e «3xl0 2 g/cm 3 . Approximating the core with a 
constant density of p He =10 2 g/cw? one gets a central He 
core of radius R He «5xl0 10 cm (0(R )). The remain- 
ing 45M Q of hydrogen will continue accreting DM (with 
a spin-independent cross section, while the helium core 
accretes it with a spin-dependent one) . In Table [T] we 
report the quantities inroduced in Section [2] as for the 
75Mq hydrogen star, as described so far, for a neutralino 
with mass m x =100GeV. The luminosity of a 75M star 
is of the order of 1O 6 L «10 39 erg s _1 throughout all 
its life, until the ox ygen burning stage, as for instance 
in (|Marigo et al 'Oil ). Our estimates show that a Dark 
Luminosity 10 37 erg s _1 < L x <10 43 erg s _1 is achieved 
for DM densities in the range [10 9 GeV cm" 3 , 10 15 GeV 
cm ]. These values are extremely interesting, as they 



comparable (actually exceed, at the cross section 



-upoer limit) the stellar luminosity during the Main Se- 
-qw 

— Transport effects: an upper limit can be obtained 
by using Eq. [5j the highest temperatures WIMPs can 
achieve correspondence to the escape velocity where most 
of the WIMPs are concentrated, namely within the ra- 
dius r x . By setting the WIMP temperature with the 
escape velocity at r x , and number densities n x and n p 
for WIMPs and protons respectively, one gets values fo 
e x as reported in the Tables. The DM transport effects in 
the core seems to be negligible for the whole neutralino 
density and mass range studied: the stellar luminosity 
L* = 10 39 erg/s needs an efficiency e* «10 9 erg s _1 cm -3 , 
if one assumes a nuclear core of 10 10 cm, two orders of 
magnitude more than upper limit value of e x achieved 
with the highest DM density considered in our range. 
During the helium burning, t x is to be read in a very 
indicative way: the higher DM density accumulated dur- 
ing the MS may simply consume itself until it reaches 
a lower concentration at the equilibrium in each part of 
the star. However, we find it interesting that a capture 
rate equivalent to the helium core is achieved by the huge 
hydrogen envelope surrounding it . 

In Table [2] we show results for low-mass stars; al- 
though the formation of low-mass stars within the Pop- 
ulation III it is not the currently favored hypothesis, 
state of the art simulations cannot totally rule out a 
low-ma ss peak in the IMF. Som e authors, as for in- 
stance (|Choudhurv fe Ferrara '06T ). do actually invoke a 
low-mass, Salpeter-like, Populationlll. 

These values are very interesting, too: as already men- 
tioned, other authors have already studied the evolution- 
ary behavior of low-mass stars in presence of a DM an- 
nihilation source, and more detailed calculations are ex- 
pected to be available soon. We only stress that for low- 
mass stars the behavior of metal-free objects is not so 
dramatically different as for massive stars, and that re- 
sults obtained for metal-rich, "usual" , stars can be used 
at least for a first approximation understanding of the 
behavior of low-mass Population III stars in presence of 
a relevant Dark Luminosity. 

4. DISCUSSION 

Without entering detailed calculations, we wish to fo- 
cus the attention on several issues that could be raised 
by an "unorthodox" behavior of massive PopIII stars, 
hoping this Letter will be used as a starting point 
for more detailed and extensive analysis. The forma- 
tion of primordial stars has to be carefully investigated 
in the presence of Dark Matter annihilation, accord- 
ing to the preliminary, interesting results obtained by 
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(Spol var et al 'Oct ); if slowed down, the cloud might start 
capturing DM at an early stage, and Dark Luminosity 
might start playing a role even before the formation of 
the actual star, if the cloud collapse is not stopped by 
this process. Assuming that a primordial star forms 
in a "standard" fashion, and lately develops Dark Lu- 
minosities of the order of our upper limits, there is no 
telling, without actual simulations and detailed calcula- 
tions, what the behavior of the object could be (we recall 
that for the 75M Q our fiducial value for X>=l(T 32 GeV 
s/cm 2 is L x =10 40 erg/s, one order of magnitude greater 
than the stellar luminosity). In the likely case that the 
actual value of uq is not its quoted upper limit, still some 
interesting scenarios can be envisioned: PopIII stars are 
thought to explode as PISNe, due to the high entropy 
and temperature of the core at the time of oxygen burn- 
ing; what if a Dark Luminosity comparable to the stellar 
one can reduce both entropy and temperature of the star, 
thus leading it through the oxygen burning phases with- 
out dramatic instabilities? What would the effects of 
an L x only comparable with the required stellar lumi- 
nosity over the short pp phase at the beginning of the 
hydrogen burning forced by the lack of CNO elements, 
given the very slow pp energy production sensitivity from 
temperature? At the same time, stellar evolution in pres- 
ence of an additional energy source might lead to very 
different nucleosynthetic yields; the odd-Z element defi- 
ciency has been long considered one of the "signatures" 
of PopIII stars, what if different temperatures (plus an 
additional electron source due to neutralino decay) would 
result in a different peculiar signature? What would hap- 
pen to a star which is supported mainly by DM burning, 
and would survive long enough to follow its host halo 
in a merger, until DM environmental conditions change? 
PopIII stars are thought to be the first engine for the 
Reionization of the Universe, as they are a huge source 
of ionizing photons as a consequence of their very high 
surface temperature; can Dark Luminosity modify it, and 
thus result in a modification of the Reionization models? 

It is straightforward to ask whether a clear signa- 
ture of this process could be recognized. High energy 
neutrinos produced by DM annihilation inside the Sun 
are expected to be observed in the future by IceCube, 
(jMena et al. '07l ): the same process would lead to a 
diffuse high energy neutrino background from the first 
stars. However, a massive star burning DM for t* wlO 6 
years with L x = 10 40 erg/s, would release a total energy 
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E x «10 54 erg; let us take as an upper limit that all of 
this energy ends up in neutrinos (the only annihilation 
ch annel being — > vis), adopt the formalism described 
in (|Iocco et al. '08f ) and th e same stellar formation ra te 
(the "fiducial" one as in (jChoudhury fc Ferrara '05l )). 
and assume the formation of a single stellar object per 
halo. One gets a flux at the Earth, at the peak en- 
ergy E v =m x /(l+z)<10GeV (m x = 100GeV) of $ <1CT 5 
GeV cm -2 sr _1 s _1 (z >10 being the "central" redshift 
of PopIII episode); this flux is buried several orders of 
magnitude below the diffuse atmospher ic neutrino back - 
ground at the same energy as from (jEvoli et al. '07l ). 
thus making neutrinos a non efficient tool for studying 
this process. 

5. PRELIMINARY CONCLUSIONS 

Within a WIMP dominated DM scenario, with struc- 
ture formation taking place in a "standard" scenario, 
primordial stars accrete DM much more efficiently than 
most of modern, galactic stars, mainly due to the pe- 
culiar conditions of the halos where they form. Within 
a large region of the relevant parameter space (DM ve- 
locity dispersion, density and WIMP-baryon scattering 
cross-section) the energy deposited inside the star by the 
accreted DM annihilation is comparable to (and even ex- 
ceeds) the stellar, nuclear luminosity. This raises several 
questions about the real nature of primordial stars, and 
about whether their behavior as stellar objects is dictated 
by baryons only. 

Whereas describing the effects of DM "burning" on the 
first stars would, at this stage, be pure speculation, we 
aim with this Letter to stimulate discussion and interest 
on this topic, and qualitatively propose scenarios whose 
validity will have to be checked in the future. How- 
ever, we think the questions raised by our preliminary 
estimates are extremely relevant for our understanding 
of first object formation and evolution in the Universe, 
which is becoming a much more complicated puzzle than 
previously expected. 
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